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Four dimensional gravity in the low energy limit of a higher dimensional theory has been expected 
to be a (generalized) Brans-Dicke theory. A subtle point in brane world scenarios is that the system 
of four dimensional effective gravitational equations is not closed due to bulk gravitational waves 
and bulk scalars. Nonetheless, weak gravity on the brane can be analyzed completely. We revisit the 
theory of weak brane gravity using gauge-invariant gravitational and scalar perturbations around 
a background warped geometry with a bulk scalar between two flat branes. We obtain a simple 
condition for the radion stabilization in terms of the scalar field potentials. We show that for 
general potentials of the scalar field which provides radion stabilization and a general conformal 
transformation to a frame in which matter on the branes are minimally coupled to the metric, 4- 
dimensional Einstein gravity, not BD gravity, is restored at low energies on either brane. In contrast, 
in RS brane world scenario without a bulk scalar, low energy gravity is BD one. We conjecture that 
in general brane world scenarios with more than one scalar field, one will again encounter the 
situation that low energy gravity is not described by the Einstein theory. Equipped with the weak 
gravity results, we discuss the properties of 4d brane gravitational equations, in particular, the value 
and sign of 4d Newton’s gravitational coupling. 


I. INTRODUCTION 

In the recent development of string/M theory [|p, branes have been playing many important roles. The idea that 
our universe is a brane in a higher dimensional spacetime has been attracting a great deal of interest. 

Although the idea of the brane world had arisen at a phenomenological level already in 1983 |Q, it is perhaps 
the discovery of the duality between M-theory and Eg, x Eg heterotic superstring theory by Horava and Witten ||] 
that made it more attractive. It actually gives the brane world idea a theoretical background: by compactifying six 
dimensions in the 11-dimensional theory, our 4-dimensional universe may be realized as a hypersurface in 5-dimension 
at one of the fixed points of a S\IZ 2 . The 5-dimensional effective theory after compactification of the six dimensions 
by a Calabi-Yau manifold can be obtained, e.g. Q. 

Randall and Sundrum proposed two similar but distinct phenomenological brane world scenarios |^,^. In both 
scenarios the 5-dimensional spacetime is compactified on SxjZ^ and all matter fields are assumed to be confined on 
branes at fixed points of the SxjZ^ so that the bulk, or the spacetime region between two fixed points, is described 
by pure Einstein gravity with a negative cosmological constant. The first scenario Q deals with a three-brane with 
negative tension as our universe. It was shown that the hierarchy problem may be solved in this scenario by a large 
redshift (warp) factor, between our brane and another, hidden brane. In the second scenario a three-brane with 
positive tension is considered as our universe and four-dimensional Newtonian gravity can be realized on the brane 
in the limit of infinite orbifold radius. 

The distance between two branes (radion) in the phenomenological RS scenario, as well as, more generally, moduli 
fields in superstring theory, should be stabilized. For this purpose a bulk scalar field 41 with certain potentials in the 
bulk and on the branes can be introduced The choice of the bulk 1/(4') and brane potentials V±('I') must be 
consistent with the 5d warp geometry . In the context of the supergravity realization of the RS scenario one has 
to consider bulk scalars of the theory with an effective scalar field potential |M,0. 

Thus, there is a broad class of brane world scenarios with a 5d bulk between two end-of-the-world branes and a 
bulk scalar field with bulk and brane potentials. This is a basis for a brane cosmology construction in which our 4d 
expanding universe is located on the visible brane. 

It is expected on the general ground that dimensional reduction of a multi-dimensional theory in the 4d low energy 
limit yields a (generalized) Brans-Dicke (BD) theory. Often 4d effective theory is derived by the integration of the 
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action with respect to the extra dimension, f dw, as it customary in Kaluza-Klein compactification. 

In the brane world scenario one can derive effective 4d dimensional Einstein equations on the brane, starting with 
the full 5d bulk equations. Although at the qualitative level the expectation that you will obtain an effective BD 
theory is justified, many quantitative details are lost. The integration f dw must be treated with special caution in 
cosmological models with brane collisions |13|. 

Widespread expectation is that BD scalar is a manifestation of a massless bulk scalar degree of freedom. If moduli 
are stabilized and no massless scalars left in the bulk, an expectation would be that BD theory is reduced to the pure 
Einstein gravity on the brane. Let us inspect a resent progress in respect to these ideas. 


II. EFFECTIVE EINSTEIN EQUATIONS, WEAK GRAVITY AND 4D GRAVITATIONAL COUPLING 


One of the most interesting results of the brane world cosmolo gy is the derivation of the 4d effective Einstein 
equations at the brane, in the scenarios with 00 or without a bulk scalar field. In the models without a 

bulk scalar the 4d effective Einstein equations on the brane have the structure 0] 

Gflu = 4 “ ^'^G , ( 1 ) 

where comes from the brane matter, A 4 is the 4d cosmological constant defined by A 4 = — 3/^q, the 

correction ir^i, include a bilinear combination of Tp,y, and Ep^ is the projection of the bulk Weyl tensor. Here, A and 
Iq are the brane tension and the bulk curvature scale (AdS radius). The coupling Gn which is interpreted in as 
the 4d effective Newton’s gravitational constant is given by 


A8n 


( 2 ) 


Thus gravity on a positive tension brane is described by 4-dimensional Einstein gravity at low energy provided the 
effects of bulk gravitational waves can be neglected. On the other hand, Gn is negative at the negative tension brane. 
In the models with a bulk scalar field (|^) is generalized to more complicated equations 00 


Gpi, — —A4{A>)gpi, + 87rG7v('I')T^i/ -I- K^iTpi, — Ep^, H— ^Tp^{’A>) , 


( 3 ) 


where Tp,y{A>) is the brane contribution of the scalar field and the effective 4-dimensional cosmological constant A 4 
now depends on the scalar field dt in a non-trivial way. Equations (i are not closed and must be supplemented by 
the bulk equations for the scalar field with proper boundary conditions. What is interpreted in [ 00 as 4d effective 
Newton’s gravitational constant is now expressed by the scalar field potential 14(4') on the brane as 

GnW = , (4) 

which is generalization of (^). One might jump to the conclusion that there is anti-gravity on the brane with negative 
potential U < 0 0 ]. 

To address the issues of the form of the low-energy theory and the value and sign of Gat it is instructive to consider 
weak gravity on the brane in a brane world scenario with fiat brane but warped 5d bulk geometry with scalar fields. 
We expect 4d Einstein-BD equations for the linearized metric perturbation hpi, around a fiat background r]pi, at the 
brane 


(□ -I- AAT) hpiy = — IGttGat 



B 


-Tv, 


'/J.2A 


2B 


(^dpd,T - 



( 5 ) 


where hpi, is defined by hpi, = hpu — hrjpv/2 with the gauge condition d'^hpi, = 0. We have introduced a factor B 
to describe the Brans-Dicke type correction to the Einstein tensor gravity. The usual BD parameter to is introduced 
in the combination B = 2 ( 3 + 2 ^) • high-energy correction. Actually, it can be shown that the high energy 

correction AAT is consistent with the so called higher derivative gravity, AAT ^ 0|. In low energy approximation 

we neglect the high-energy correction, AAT 0. 

Weak gravity in the Randall-Sundrum brane world (without bulk scalar) was investigated by Garriga and Tanaka 0 ] 
and they showed that linearized Brans-Dicke theory is realized on the brane with 
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G% = 


3±(i/Zo 


WttIo sinh(d/Zo) 


B± = , 


( 6 ) 


where d is the inter-brane distance. Einstein theory is restored only in the limit of an infinite orbifold radius d, to on 
the positive tension brane becomes infinite, 5+ = 0. On the other hand, formula (|l]) of Shiromizu, et.al jl^] shows 
that gravity on a positive tension brane is described by 4-dimensional Einstein gravity at low energies provided the 
effects of bulk gravitational waves can be neglected. 

In order to understand the two different-looking results (|^) and (|^) in a unified manner, let us briefly review 
linearized Einstein and Brans-Dicke theories around a 4-dimensional Minkowski background. Eirst, because of the 
Poincare symmetry of the background it is possible to Eourier expand all perturbations and to classify them by 
irreducible representations of the little group, see Appendix for the details. There are three relevant classes: scalar, 
vector, and tensor perturbations. Next, we can choose a suitable gauge or construct a set of gauge-invariant variables. 
Third, we can write down linearized equations for each class separately. It is easily shown that vector perturbations 
should vanish. The characteristic feature of Einstein theory is that the coupling constant betwee n gra vity and matter 
is the same for scalar and tensor perturbations and given by Newton’s constant (see subsection V A). On the other 
hand, in Brans-Dicke theory the universality is spoiled by the existence of a scalar field included in the theory since 
the scalar field couples with the scalar perturbations only. The Brans-Dicke parameter determines how the scalar field 
couples to scaler perturbations. Hence the gravitational coupling constant for scalar perturbations is a combination 
of Newton’s constant and the Brans-Dicke parameter while that for tensor perturbations is Newton’s constant itself. 
To see this, let us set AK = 0 in (||) to obtain the linearized equation for 4-dimensional BD theory and split 

-I- — hrif^^/A, with the transversal-traceless 

LL 

fll' 


into tensor (transverse-traceless) and scalar (trace) parts, 


conditions upon /iJJ. 


Here, is given by 


^ (df,d^h - 


(7) 


Accordingly, equation (||) is split into two parts, tensor and scalar ones 


TT 


= -leTrGjvT),,. 


Dh = I67rGjv(l - 2B)T . 


( 8 ) 


Let us apply this to the issue of weak gravity on branes in the Randall-Sundrum scenario. We can easily show by 
using the 5-dimensional Einstein equation that there are no scalar-type perturbations in the bulk except for a zero 
mode called radion. The zero mode does not affect the value of the 4-dimensional gravitational coupling constant 
for scalar perturbations G^(I — 2B±) in (H) since D/i = 0 for a zero mode by definition and T = 0 for a zero mode 
by the conservation equation. Hence, we can pretend as if there were no scalar-type perturbations in the bulk, and 
the results (0), i) of Shiromizu, et.al ||l5| without ambiguities determines the 4-dimensional gravitational coupling 

constant for scalar perturbations Gn = G^(I — 2B±) = ±g^. This follows from (^) and from (||) where one has to 
use flat brane tuning A = 

On the other hand, the coupling constant for tensor perturbations cannot be determined directly since there are 
tensor-type perturbations in the bulk manifested in In other words, the system of equations ( 0 ) is not closed 
in respect tensor modes, and shall be supplemented by information about bulk Weyl tensor. Actually, the existence 
of tensor-type perturbations in the bulk E^i, shifts the 4-dimensional gravitational coupling constant for tensor 
perturbations from Gn as it defined in ( 0 ), (^), see also The gravitational coupling for tensor mode is given 
by Gat from (^. As a result, we have different 4-dimensional gravitational coupling constants for scalar and tensor 
perturbations and, thus, the linearized Brans-Dicke theory. 

Let us discuss more general case of the brane world scenario with bulk scalar fields. 

We begin with the case of a single bulk scalar field dt. Equation (^) describes effective 4d gravity on the brane. 
We just learned that bulk gravity is important to reproduce correct weak gravity at the brane. Now we have to 
understand the effect of bulk scalar contribution T^iy(4') and A4(4') to the 4d gravity in equation (||). 

In the analysis of weak gravity one of the essential differences between models with and without a bulk scalar field 
is that there are non-trivial scalar-type perturbations in the bulk for models with a bulk scalar field, whereas these 
perturbations vanish for models without a bulk scalar field. Hence, the existence of the scalar field can shift the 
4-dimensional gravitational coupling constant for scalar perturbations from what is given in (^. On the other hand, 
we already know that the existence of tensor-type perturbations in the bulk shifts the 4-dimensional gravitational 
coupling constant for tensor perturbations. Hence, our question now is whether these shifts in the gravitational 
couplings for scalar and tensor modes due to bulk gravity and bulk scalar T)iiy(4') are the same or not. If they 
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are the same then the effective gravitational theory on the brane is pure Einstein theory (with a shifted Newton’s 
constant, which differs from (ft)). If they are not the same then the effective theory is still Brans-Dicke theory, but 
still shifted from (|). 

The issue of weak gravity in the generalized braneworld scenario with a single scalar field was already investigated 
by Tanaka and Montes They obtained very intriguing result that 4-dimensional Einstein gravity is restored at low 
energy. In other words, according to the above arguments, their claim is that shifts in the 4-dimensional gravitational 
constants for scalar and tensor perturbations bring them to the same value. 

The analysis of [0 can also be generalized to include non-minimal coupling (which we denote as e “±('*'±1) of the 
brane matter to the metric in the 5-dimensional Einstein frame (to reflect the spirit of dilaton-matter coupling). This 
non-minimal coupling was included in the 4d Einstein equations in ||l^ , |l^ , where it was assumed that matter on the 
brane is minimally coupled not to the metric in the 5-dimensional Einstein frame but to a conformally related metric. 
It is evident that the induced metric on the brane in the minimally coupled frame should be considered as a physical 
metric on the brane. 

Since the conformal factor depends on the bulk scalar field, metric perturbations in the minimally coupled frame 
have an extra contribution from scalar perturbations. Hence, from the above arguments about the universality of 
gravitational coupling constants for scalar and tensor perturbations in 4-dimensional Einstein gravity it is expected 
that the restoration of 4-dimensional Einstein gravity may depend on the choice of the minimally coupled frame. 

In the case of a single bulk scalar field we obtain 


G 


± _ 
N — 


IGtt J ’ 


B±=Q . 


( 9 ) 


Here, is the warp factor in the 5-dimensional Einstein frame, whereas the warp factor in the minimally coupled 

frame is The value of G% is always positive and differs from ®). We conclude that the coupling with the 

brane energy-momentum tensor alone in the effective Einstein equation (^ cannot be interpreted as the 4d Newton’s 
constant. The correct value of the Newton’s constant is given by (D- In particular, there is no antigravity at the 
negative tension brane. 

The purpose of this paper, on technical side, is to perform a complete analysis of weak gravity in the generalized 
braneworld world scenario in which there is a bulk scalar field. We take into account the fact that matter on the 
brane is not minimally coupled to the metric in the 5-dimensional Einstein frame in general and allow an arbitrary 
conformal transformation to the minimally coupled frame. Despite the above expectation, our conclusion is still the 
same as that of Tanaka and Montes: 4-dimensional Einstein gravity is restored on branes at low energy. We formulate 
the necessary and sufficient condition for brane gravity at low energy to be the Einstein theory. We also will try to to 
understand the connection between the bulk scalar degrees of freedom and the character of the brane gravity. Indeed 
Einstein, not BD, gravity at the brane is restored in the weak gravity for a dilaton potentials in the bulk and at 
stabilized branes. Moreover, this result can be extended for the second order of perturbative gravity within a certain 
class of symmetric perturbations |2|] . 

It is conjectured that for two or more bulk scalar fields (unless their potential has special properties), the weak 
gravity at the brane may be again BD gravity. We can summarize that for models without a bulk scalar field the 
parameter B in the low-energy gravity theory (||) is H . In the case of a single bulk scalar field B = 0. For 

n > 1 scalars in the bulk we may again expect a nonzero B. It is intriguing that the BD parameter, which can be in 
principle constrained by observations, depend on the details of the bulk theory. 

How one can understand this result? Suppose there are n scalars in the bulk. In the weak gravity limit we still 
expect the effective gravitational theory of the form (Q). In other words, we expect the BD theory with a single 
effective BD scalar on the brane. Thus, the effective BD scalar on the brane, in some sense, is a composite object 
due to superposition of contributions from the bulk scalars. For the case without scalars in the bulk, but the radion, 
we have the non-zero value (D of B. For the case of a single bulk scalar the superposition of its contribution and the 
radion contribution to the brane scalar gravity are precisely canceled so that B = 0. In the case of more than one 
bulk scalars the precise cancellation is not always expected. _ 

The rest of the paper is rather technical and is organized as follows. In Sec. |h| we describe our model of the brane 
world and define the physical energy momentum tensor of matter fields on the branes. In the section we also specify 
the background around which we analyze perturbations in Sec. IV. In Sec. ^ we compare gravity in the brane world 
with that in 4-dimensional Einstein gravity. We show that 4-dimensional Einstein gravity is restored at low energies 
on both branes provided there is no tachyonic mode with a boundary condition corresponding to vanishing stress 
energy of matter fields on the branes. Finally, Sec. VI is devoted to a summary of the paper and discussion. 
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III. MODEL DESCRIPTION AND BACKGROUND 


In this section we describe our brane world model and specify a background. In the next section we investigate 
perturbations around this background to analyze weak gravity in the brane world model. 



FIG. 1. Sketch of the brane world geometry. 


Let us consider a 5-dimensional spacetime j\4 with topology AI 4 0 IZ 2 , where AI 4 represents 4-dimensional 
spacetime. It is easily understood that two timelike hypersurfaces corresponding to two fixed points of jZ^ 
can be sources of gravity. In fact, because of the Z 2 symmetry, the first derivative of the fields can be discontinuous 
across each of the fixed points and thus, the second derivative can include a delta-function. Since Einstein equation 
is a second order differential equation, inclusion of the delta-function in the second derivative of the fields implies the 
existence of gravitational source at each of fixed points. Since we would like to consider E± as a world-volume of 
a brane and would like to think of one of branes as our 4-dimensional universe, we assume that the 4-dimensional 
intrinsic geometry on E± is regular. On the other hand, the 5-dimensional geometry is not necessarily regular on E±. 
The 5-dimensional region off both of E± is usually called the bulk spacetime and we denote it by M.h- Hence, 

M=M6UE+UE_. (10) 

To describe the hypersurfaces E± we use the parametric equations 

E±: x"^ = Zf( 2 /±), (11) 

where {M = 0, • • •, 4) are 5-dimensional coordinates in A4 and each y± denotes four parameters {y^} (/r = 0, • • •, 3) 
playing the role of the intrinsic coordinate system on E±. See Figure ^for illustration. 

We consider a dilaton type theory with a single scalar field to describe gravity and assume that all other matter 
fields are confined to the hypersurfaces E±. Hence, in the following arguments we assume that the system is described 
by the action [| 


Itot — Ieh -\- Ii 


matter 5 


( 12 ) 


where Ieh is the Einstein-Hilbert action 


^For simplicity we do not consider the boundary of M 4 , but it is easy to take it into account by imposing suitable boundary 
conditions and introducing boundary terms appropriate for these boundary conditions. 
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( 13 ) 


IeH = TT^ [ d^Xy/^R, 

2^5 Jm 


Iq, is the action of a scalar field 'I' 

= — d^xy/—g 

dMb 




“ X! / d'^ya^J-q^Vai^a) 


-L -'So 


G—±. 


and Imatter IS the action of matter fields living on the branes 

^matter — / d^2/o-^±[9±A»i/,matter], 

a=±-'s- 


(14) 


(15) 


Here, 4*^ and q±^v represent the pullback of the scalar held 4' and the induced metric on and q±fiu is the physical 
metric on T,± specihed below. The physical metric q±fj,i, is actually ‘physicab in the sense that because of the form 
( p^ ) of the matter action, all matter helds living in the 4-dimensional manifolds E± feel the geometry described by it. 
The induced metric q±fiu and the pullback 4'j- of 4^ on S± are dehned by 

Q±fJ.viy±) ) 



^±(2/±) '^\x=Z±{y±) ’ 

(16) 

respectively, where 

M dZ^ 
dy^ ■ 

(17) 

We assume that the physical metric q±fj,i, 

is conformally equivalent to the induced metric q±^,v'- 



q±y.v = exp[-a±(4'±)]g±^^. 

(18) 


where a± is a function of 4^±, respectively. 

As shown in ref. the action above gives a consistent variational principle including variations of the positions 

of the hypersurfaces as well as the metric gMN, the scalar held 4>, and other matter helds. Specializing to our 
model with the Z 2 -symmetry, the equations of motion are as follows. 


Gmn = kITmn, 

2K±,. = , 


and 


4- 


M 

;M 


- 17'(4') = 0, 

2a±4'± = H;(4'±) + ia'±(4'±),S±, 


(19) 


( 20 ) 


where the semicolon denotes the covariant derivative compatible with the metric gMN, Gmn is the Einstein tensor 
for the metric gMN, Tmn is the stress energy tensor of the scalar held in the bulk, and prime is a derivative with 
respect to 4'. 


Tmn = Om'^Qn'^ — gMN 


]^g^'^'dM''i’dN''i> + U{^) , 


( 21 ) 


and 3j_4'± are, respectively, the extrinsic curvature of and the normal derivative of 4^ on dehned by 

^ N[ N n I 

R±ii.v — 2 , 

5x4-±= (22) 

is the unit normal to the hypersurface S± directed towards the bulk S±yv is the surface stress energy tensor 
given by 
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( 23 ) 


= Q^fipQ^utrS^ , 

Sf^"'=—^=jJ -— [ d'^y±C±[q±fj,^,mutter], 

V-d± Sq±pu Jt.± 


and 


S± = = -4y±(^±) + e-2“±('^±)5±, 

5± = r''S±^,.. (24) 


Note that can be called the physical surface energy momentum tensor since it is defined in terms of the matter 

action and the physical metric only. Hence, when we seek an effective 4-dimensional gravitational theory on the 
branes, we have to write the equations in terms of S±fj,i, instead of S±fiu- 

Now let us specify the background {qmn = , T±^i, = T®,,) around which we shall 

investigate perturbations in the following sections. We consider a background with 4-dimensional Poincare symmetry: 


dx^ = e ^"^^^\r]^,^dx^dx‘' + duP'), 

^(0) ^ ^(0)(y;)^ 

y(e)p _ U 

- y±^ 

r7{0)w 

=w±, 

c(0) _ Q 


(25) 


where x^ (p, = 0, • • •, 3) represent the first four of the 5-dimensional coordinates x^ {M = 0, • • •, 4) in A4, w represents 
the fifth coordinate x^, and w± are constants. Here, we have redefined V± and L± so that vanishes Note 
that when we consider perturbations around the background as in the next sections, x^ on E± does not necessarily 
coincide with y^. 

The equations of motion for the background are as follows. 

3i-b3i^ = 

3H - gi^ = 2«;^e"^'^f7(«'(°)) (26) 


and 


s^\i/(o) 


w—w± 


w—w± 



(27) 


where dots denote derivative with respect to w and = 4'^°)(wi). Note that the sign in the right hand side of 
( p^ ) is due to the sign in the following expression of the unperturbed unit normal to S® directed towards Mb- 

n'l^^dM = (28) 

Here, we assumed that W- < w+ and that the bulk is the region w- < w < w+. The induced metric and the physical 
metric on E± are 


„(o) 

JO) 


_ „- 2 A(w±) 


( 29 ) 


^This is possible because of the Poincare symmetry. 
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IV. PERTURBATIONS 


In this section we investigate perturbations around the background specified in the previous section. Namely, we 
consider 


9m N = g^MN + ^9MN, 

= ^f(o) + S'if, 

+ sz^, 

= + (30) 

The unperturbed functions specify unperturbed hypersurfaces and the perturbed functions {Z^} specify 

perturbed hypersurfaces E±. As pointed out in ref. [psf there are two kinds of gauge transformations, one in the bulk 
and another on the branes. The former is called 5-gauge transformation and is of the form 

X^^x^+^^ix). (31) 


The latter is called 4-gauge transformation and is of the form 

2/^^y^+Cl(j/±). (32) 

Following ref. [^, it is straightforward to calculate perturbations of the induced metric q±^i, and the extrinsic 
curvature from 6gMN and Z^. The result for a general background is 


where 


g±fj,u — g±^v T ^g±p.i'i 


Sq±t,u = {6gMN + Csz^gM^)^ 




( 0 )M( 0 )N 


n''^>^\SgMN + 2SZ±M-N)K^±l 


(0)L (0)M (0)iV 


“^STlmn + SZ±l;mn + SZ±l-nm + (Rl' 


(0) , n(0) 

L'MLN ^L'NLM 


L' 


)SZ^ 


(33) 


(34) 


Here, C denotes the Lie derivative defined in the five-dimensional spacetime, the semicolon denotes the covariant 
derivative compatible with the background metric 5^^, R^l^mln i® Riemann tensor of g^^pf and ST^mn = 
(ll2){5gLM-,N + 5gLN-,M — ^gMN-z)- The right hand sides of the above expressions for perturbations Sq^^, and SK^i, 
can be evaluated on the unperturbed hypersurface S® at cc = Z^\y±). It was shown in ref. that 6q±fj,i, and 
5K±^i, are invariant under the 5-gauge transformation and transform under the 4-gauge transformation as 


^ (35) 


where C denotes the Lie derivative defined in the 4-dimensional manifold Trivially, SS±fj_i, is also invariant under 
the 5-gauge transformation and transforms under the 4-gauge transformation as 

^ (36) 

Applying these formulas to our background, we obtain the following expressions for 5q±^i, and 5K±^i, 


^q±HL' 

5K±^u 


4 “ 2 ^ ^^^gwu di/Sg^ufi ^g^u^gww 

-2A{d^5Z^, + d^5Z^,) - 2d^dJZ±^ + 2e^ie-^y g^JZ±, 


( 37 ) 






where 5Z±m = 9m\^^±- 

Similarly, we obtain the following expressions for perturbations of 'f'i and 

+ 6^±, 

(38) 

where 


= S'l> + e‘^^¥°^SZ. 




(59^'hj- = =p 


+ 2e^(5^ + SZ±^ 


(39) 


Again, the right hand side can be evaluated on the unperturbed hypersurface S® at a; = z‘f\y±). It is easy to 
calculate perturbations of the physical metric q±fiiy on by using the following relation 


( 0 ) , 




V, (0)^-aWg-2A± 




±) 


(40) 


where A±, a® and represent A(w±), q;±( 4>®) and a(j_('I'£^), respectively. 

We have the following 4-gauge transformation. 


^ ^uCizfn 

5K±^i, (5(invariant), 

5S±^i, (55±^;.(invariant), (41) 

where C±m = q±lX±, 

5k±^, = + K^il^5q±p^), (42) 

and To see the 5-gauge invariance of note that we have redefined V± and C± so that 

^±l.u vanishes (see eq.(p5|) and the footnote after that.). The perturbed junction condition can be written in terms 
of these variables as 

SS±f,, = + e-“±^5±^., (43) 


2SK±., = -kI 6S±., - 


and SS± = q±'^'^''SS±^i,, where represents ^(.('I'^y The scalar field matching condition is written as 

=F -h (44) 

where and 6S± represent and q±^'^‘'SS±fj,v, respectively. 

On the other hand, the 5-gauge transformation is 

^ 

Sg^w Sgfj.1^ - - d^^w, 

^gww ^ ^gww 2e (e j 

SZ±M SZ±M + ^M, 

S^^J ^ S'i/- (45) 


where 

Since the background has 4-dimensional Poincare symmetry it is convenient to expand the perturbations in terms 
of scalar, vector and tensor harmonics in 4-dimensional Minkowski spacetime: 
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SgMNdx^dx^ = {h(T)T(T)fiv + h(LT)T(LT)fii^ + h(^LL)T(LL)^iv + h(^Y)T(Y)tJ.u)dx'^ dx‘' 

+ 2{h(^T)wV{T)fj. + d(L)wV{L)fj.)dx^dw + ^ww Ydw^, 

S^i> = ^Y, 

5 Z±Mdx^ = (z±(T)V(T)^ + z±(^L)V(L)^)dx^ + z±yjYdw, ( 46 ) 

and 


ds±^„ — T±t^T)T(T)iii^ + T±{LT)T{LT)^iv + T±{LL)T{LL)iiv + T±{Y)T{Y)iivi (47) 

where T(j’ [^j’,LL,Y)fiuj y{T,L)^ and Y are harmonics in 4-dimensional Minkowski spacetime (see appendix |A| for their 
definition), the fc-dependent Fourier coefficients and ip are functions of w only, and the 

other coefficients z±(T,L)i z±w, and t±i^t,lt,ll,y) constants. We omitted the integration with respect to k. Let 
us expand the gauge parameters and too. 

^Mdx^ = {i(T)V(T)^i+ ^(L)V(L)^j)dx^^ + ^yjYdw, 

C±M = C±(T)4^(r)/i + C±(l)4^(l)^> (48) 

where and are functions of w only and C±{t,l) are constants. 

Using the harmonic expansions (^^, we can obtain the corresponding harmonic expansion of Sq±^,j, K±^v, 
and (59i'I'±. The result is 


dq±^lv — ^±(T)T{T))iv + f^±{LT)T{LT)iii^ + <^±{LL)T(LL)iil^ + f^±(Y)T{Y)fii^, 
5 K±^ 1 , = k±t^T)T(T)fj.v + k±(LT)T(LT)ni' + k±(LL)T{LL)fj.v + k±(Y)T{Y)^v, 
i5^'± = ip±Y, 

5 d±_'i>± = ip±iY, 


( 49 ) 


where ^-dependent Fourier coefficients are 


( 0 ) 


^±(T) = e “± F'(t)j 
_ ( 0 ) 

0 '±(LT) = e “± <(>±(T), 

_ ( 0 ) 

0'±(LL) = e “± <()±(L), 

( 0 ) 

ct±(y) = e “± 

k±{T) =T-j^e.~^{e^^F{T))P 
k±(LT) = 


F - 2A(j)±^ - -ri^''kf,ky(j)±i^L) - a'±\e + (()±„4'(°)) 


k±(LL) = F-e ^±iu, 

k±(Y) = ^e-^{e^^Fy + (e^) U^^ + 

y)±=ip + 

1 


-e^r]^''kf,kt, - 2(e^)' 


(l>± 


^±y=T2 




±10 


( 50 ) 


Here, the right hand sides of 
invariant variables defined by 


are evaluated at w = w± respectively and have been written in terms of 5-gauge 


F(t) — k(T), 


Fn, = h. 


(T)v 


„-2A^2A 


(e h(LT))', 


F — h(Y) + 2AXjju + kf^ki,h(^]^p, 

F — h — 2e~^(e^X V 

ww — '^WW J 5 
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0±(T) = Z±(^T) + , 

</>±(L) = Z±{L) + h(LL)\^^^^ , 

4^±w — Z±uu + , (51) 

where It is easy to show that these variables are actually 5-gauge invariant, by using 

the following 5-gauge transformation. 


k{T) 

h{T), 


h(LT) 

i 

si 

1 

vn 

7s 


h(LL) 

1 

T 


hr) 


> 

h{T)w 

^ h{T)w - e“^"^(e^"^^(r))', 


k(^L)w 

^ h{L)w — — e , 


kww 

—j. ?7 _ V 


Xu, 

1 

T 

(52) 

The reason why the coefficients of 6 q±^i,, 

5K±^v, '5'1'i and were 

written by 5-gauge-invariant variables 


only is that they themselves are 5-gauge-invariant as shown in ref. for and 6K±^^. As for the 4-gauge 

transformation, while 5K±^^, and i9j_'I'± are 4-gauge invariant, 5q±^i, is not. Its 4-gauge transformation is 


<^±{T) cr±(r), 

o'±(Lr) ^ (^±(LT) - C±(r)) 

<^±(LL) ^ <^±(LL) - C±(L), 
o-±(v) ^ (t±{y) + 

Hence, we can construct the following doubly-gauge-invariant variables. 


/±(T) — 0-±(T) — 


( 0 ) 

— CKm 77’ 

2 ^ t. 


(T), 


- _ 1 _ _ 

/± = cr±(y) + i^V^''kf,K(T±(^LL) = e “± 




(53) 


(54) 


Note that the coefficients of 5S±^i, are doubly-gauge invariant variables by themselves. 

Since there are many coefficients in the above harmonic expansions, let us divide these into three classes. The first 
class is the scalar perturbations and consists of coefficients of Y, V( 2 ,)^, T(^LL)iiv and T^^Y)^^v■ The second is the vector 
perturbations and consists of coefficients of V(t)^ and T(^LT)fj.v The last is the tensor perturbations and consists of 
coefficients of It is not difficult to show that perturbations in different classes are decoupled from each other 

to linear order. Hence, in the following we analyze perturbations in each class separately. Decomposition into scalar, 
vector and tensor modes is commonly used in cosmology. However, usually in cosmology we use scalar, vector and 
tensor representations of the isometry group related to the symmetry of 3d space. Meanwhile, here we will use scalar, 
vector and tensor representations of the isometry group of 4d space-time. 

Before analyzing each class in detail, here we summarize some differences between our formalism and one in the 
literature [O,^. What we have to do is exactly the same in both formalisms: we have to solve perturbed Einstein 
and scalar equations in the 5-dimensional bulk with boundary conditions given by Israel’s junction condition and the 
scalar field matching condition on the branes. Of course, this system can be considered as a set of inhomogeneous 
differential equations with Z 2 symmetry since the non-trivial boundary conditions can be interpreted as source terms 
with the delta-function form Hence, mathematically there are two approaches to this problem: one is to 

solve the differential equations with the boundary conditions (or the inhomogeneous equations) directly; another is 
to use the Green’s function method to solve the inhomogeneous equations by using the complete set of homogeneous 
solutions. The formalism in the literature |p^ , ^ is based on the second approach. In this approach we need the 
complete set of homogeneous solutions in order to construct Green’s function. The set of homogeneous solutions 
consists of the so called zero modes and Kaluza-Klein modes. These homogeneous solutions satisfy the boundary 
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conditions corresponding to vanishing matter stress energy on the branes. Accordingly, the homogeneous solutions 
have discrete spectra. On the other hand, in our formalism we treat the differential equations with the boundary 
conditions (or the inhomogeneous equations) directly. By doing so, we can Fourier transform all equations with 
respect to the 4-dimensional coordinates to reduce the problem into a set of purely 1-dimensional problems. We also 
classify the Fourier components into irreducible representations (scalar, vector and tensor modes) of the little group 
of 4-dimensional Poincare symmetry. Note that is continuous in our formalism. 

Although r]^'^kfj,ki, is continuous and there is no discreteness in it, the scalar and tensor Kaluza-Klein masses 
may still have important roles in our formalism. Here, a scalar (or tensor) Kaluza-Klein mass uikk is, as in the 
literature, such a non-zero value of ■\/—r]>^‘'k^ki, that the scalar (or tensor, respectively) perturbation equation have 
a non-trivial solution with vanishing matter stress energy on the branes (homogeneous solution). It is expected that 
rj^'^k^ki, = corresponds to a resonance point in the continuous spectrum of rj^'^k^ky. Hence, there may be 

interesting phenomena near the point corresponding to the Kaluza-Klein mass, rj^'^kfj^ki, ~ However, we shall 

not see this in this paper since we shall adopt the expansion by a parameter = Prj^'^k^k^, where I is a length scale 
of the model. Further analysis near the resonance may be an interesting future work. 


A. Scalar perturbations 

For scalar perturbations, we can use the the background equation to reduce the perturbed junction condition to 


(t)±w = TK5e“^““±”f±(ii). (55) 

Here, we have assumed that k^ ^ Q since a mode with k^ = 0 preserves 4-dimensional Poincare symmetry and 
can be included in the background. The first equation is a gauge-invariant expression of the perturbed conservation 
equation. The second equation gives the gauge-invariant perturbation (f>±w of the position of each brane in terms of 
matter perturbations on the brane. The perturbed matching condition of the scalar field at w = w±, respectively, is 


T 


„3A 




■ 2e"^(f ■ 








(56) 


where, as we already stated, and are abbreviations of K(.'('I'£^) and a'_j_('I'£^), respectively. 

For scalar perturbations, the perturbed Einstein equations in the bulk consist of the following set of equations for 
5-gauge invariant variables. 


= —2F, 


2k24'(o) ’ 


(57) 


and 


F 


2A- 


(e^y 


(e^y 


F - [2e-^(e^) - + F = 0. 


(58) 


In the following arguments, we shall consider modes with rj^'^kfj^ki, = 0. These modes represent massless fields 
on branes and are called zero modes. As we shall see in the next section, for scalar perturbations it is actually 
sufficient to analyze zero modes in order to compare the low energy behavior of gravity in the brane world scenario 
with 4-dimensional Einstein gravity. 

Setting rj^'^k^k^ = 0, it is easy to find a zero-mode solution F = F^^"^ = (e^)' of (^^. Another independent solution 
F = F^‘^'> can be easily found by using the fact that the Wronskian A = obeys the equation 

A = —HiA, where Hi is the coefficient of F in (^^. The solution to the Wronskian equation is A oc e“^^(e^) ', and 
thus 


^(2) 




dw 


A + 

-:-( 

A2 


,-3A 


Hence, a general solution is 
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( 59 ) 


F{w) = T + ^K)^ g-3AK). 

7o A{w'Y 


lo A{w'Y 

For the solution (|9|), the perturbed matching condition (|5^) of the scalar field is 


(Ci + C2B{w±) ± 2 K 2 e- 2 ^±-“™f±(ii) j i?± = 0 , 


(60) 


where 


7o ^(rc')" ^(^i') 

B± = (i + i^) ± 

= ^g-4Ap/(0)2 ^ ^yWy^iO) _ ly^(0)^„(0) 


(61) 


and A±, and {7^*'°^ are abbreviations of A(w±), a±('I'^^) and {/'(dr^^), respectively. Here we have used the 
perturbed junction condition ( ^^ and the background equations to simplify the equation. Hence, if B± ^ 0 then the 
coefficients Ci and C 2 are uniquely determined as 


Cl 

C 2 


2^1 


B{w+) — B{w-) . 

24 

B{w+) — B{w-) . 


B{w-)e 


-^^^-<'f+iLL)+B{w+)e- 


- T_ 


(LL) 


+ + t+(ll) + e - r_ 


(LL) 


(62) 


Note that there are no further conditions and that this solution represents a set of gauge invariant perturbations. 

For further analysis it is rather important that B± 7 ^ 0. The rare case of B± = 0 would correspond to the situation 
when scalar brane matter perturbations are not coupled with the bulk perturbations. It takes place, for instance 
for the RS models with constant U and V±. In general, the forms of the potentials not necessarily lead to B± = 0. 
It is interesting to check that B± = 0 in the 5d reduction of HW theory ||] with 1/(5') = exp (— 2 '\/ 2 «; 5 'I') and 
V± = ±/3-\/2exp (—•\/ 2 'v 55 '), /3 is a constant. For these potentials the radion is not stabilized in the sense that the 
inter-brane distance is arbitrary. 

Therefore, by using the formula (p^, the doubly gauge invariant variable / representing perturbations of the 
physical metric is written as 


f± 


-2A±- 




dw' 

JW — 


,- 2 A+-a™- 


t+{ll) + e 


-2A--A°'’ - 


-(LL) 


(63) 


Here, we have used the fact that B = to simplify B{w+) — B{w-)i and we have also used the background 

equation. 

Although we worked in a sector with = 0, the result (|^) can be considered as the lowest order term in the 

expansion of f± by the parameter /r = Prj^'^kfj^k^, where I is a characteristic length scale of the model. (See ( [7^ for 
the expansion of tensor perturbations.) Hence, the above result for ^ = 0 is enough to give the following expression 
of the lowest order term in the ^ expansion for a small non-zero /r. 


f± 


,- 2 A±-a™ 


Jw_ 


,-2A+-. 


+ T+(ii) 




- T^-(LL) 


0 ( 4 . 


(64) 


B. Vector perturbations 

For vector perturbations, the perturbed junction condition at w = w±, respectively, is 

B^u = TK 5 e"^±““±’r±(iT), (65) 
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and there is no matching condition for the scalar field. 

The perturbed Einstein equation in the bulk for vector perturbations is 


= 0 , 

{e-^F^y = 0 . 


Hence, the solution is 


J Ce^, for q^''ki^kiy = 0 
\ 0, for q^^k^k^ ^ 0 


( 66 ) 


(67) 


where C is a constant. For q^'^k^ki, = 0, the junction condition overdetermines the coefficient C: 


C = -4e- 


dO) . 


2-2A--a 


( 0 ) 


+ T+(LT) — 1^5^ ““ “ ''■-(LT)- 

Hence, for q^'^k^ki, = 0 we have the following relation between r+(LT) and 

-2A+-a‘'°^ - I r, 

e + +e - T-{lt) = 0. 


( 68 ) 


(69) 


As we shall see in the next section, this relation is important for restoring 4-dimensional Einstein gravity at low 
energies. For q^'^kf^k^, 0, the junction condition becomes 


■^-KLT) — T-iLT) — 0. 


(70) 


Note that for vector perturbations there is no doubly gauge invariant perturbation of physical metric on the branes. 


C. Tensor perturbations 

For tensor perturbations, the perturbed junction condition at w = w±, respectively, is reduced to 

(e^^F(r))' =±K5e^±““± t±(t) (71) 

and there is no matching condition for the scalar field. 

The perturbed Einstein equation for tensor perturbations in the bulk is 

gA = 0. (72) 

It can be shown that q^y^ki, < 0 for any solutions of the equation (|^) insofar as the boundary condition (0) 
with = 0 is satisfied i^. The proof is easy: we can show that 

q^'^k^ki, / dwe^F'^rp^ = — dwe~^^ ] < 0. (73) 

J W— J W— 

In the following arguments, we shall first consider zero modes, or modes with q^'^k^k^ = 0. Contrary to the scalar 
perturbations, we will see in the next section that for tensor perturbations it is not sufficient to analyze zero modes in 
order to compare the low energy behavior of gravity in the brane world scenario with 4-dimensional Einstein gravity. 
Hence, after analyzing zero modes, we shall investigate non-zero modes, or so called Kaluza-Klein modes by perturbing 
with respect to q^'^k^k^ around the zero modes. 

For q^’^kpki, = 0, the solution of the bulk Einstein equation ( |7^ ) is 

pW 

F^t){w) = +D2e-^^^^\ (74) 

Jo 

where Di and D 2 are constants. By using the formula (|5^), the doubly gauge invariant variable /±(t) representing 
perturbations of the physical metric can be written as 

/±(r) = -f £>2e"2^±"“±\ (75) 

^0 
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The junction condition overdetermines the coefficient Di\ 


D\ = K,ie 


„( 0 ) 


T+(T) 


2-2A--a 


( 0 ) 


-(T). 


Hence, for = 0 we have the following relation between Tj^(t) and 


( 0 ) 

r+(T) + e 


f_(r) = 0 . 


(76) 


(77) 


As we shall see in the next section, this relation is important for restoring 4-dimensional Einstein gravity at low 
energies. 

For 7 ^ 0 it is not trivial to solve ( [^ without specifying the model potentials. Hence, we shall perform 

perturbation with respect to fj. = around the zero mode solution, where I is a characteristic length scale of 

the model By expanding and f(y) as 

E(t) = + O(p^), 

r±(T) = r±(r) + + 0{n^), (78) 

we shall solve the Einstein equation ( [^ in the bulk and the junction condition 0 ) order by order. The zeroth order 
analysis is exactly the same as the previous analysis of the zero modes. Hence, the zeroth order Einstein equation in 
the bulk has the solution 


pw 

(w) = j + £>26 


-2A{w) 


where Di and D 2 are constants, and the zeroth order junction condition overdetermines the coefficient Di. 


Di = nje “+''^+(T) = 


2-2A.-a^y-[0] 




and thus we have the following relation between and . 


- 2 A+-a™-[ 0 ] 

^ ^+(T) ^ ® 


+ (T) '-(T) 

- 2 A_-a™=[ 0 ] 




Note that the zeroth order analysis leaves the coefficient D 2 undetermined. 
To first order in /j,, the Einstein equation (^ becomes 




g-3A(g2A^[l]). 


_ ^[ 0 ] _ Q 


Hence, we obtain 


^2g-3A(g2Aj^m^). ^ + D 2 B 2 (w) + D 3 , 


( 79 ) 

(80) 

(81) 

(82) 

(83) 


where D 3 is a constant and 


Bi(w) 

B2{w) 


pW pH 

/ / 

Jo Jo 

pw 


dw e 


// 3A{w") 


(84) 


By integrating ( | 8 ^ ) once more we obtain a solution of the bulk Einstein equation to this order. For the solution, the 
junction condition ( 0 ) becomes 


®The characteristic length I is determined by the background solution and determines the energy scale below which the 
low energy description is valid. 
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(85) 


DMw±) + D2B2 {w±) + ^3 = rgy). 


Hence, the constants D 2 and £>3 are uniquely determined by In particular, 

Bi{w+) - Bi{w-) 


^2l2 


D 2 = 


B2{w+) - B2{W-) . 


-2A+-a™-[l] 

e + + + e 


±(T)- 


B2{w+) - B2{W-) 


Di. 


( 86 ) 


Therefore, for small |/i| = by using the formula (p4[), the doubly gauge invariant variable /±(t) repre¬ 

senting perturbations of the physical metric can be given by 


JW— Jw+ 


„2A±+a™ 


/±(r) = -Di 


W± 


Jw_ 


-h- 


4P 


f’”’*' dw '- 

J W — 


+e- 
HT) ^ 


-2A-- 


( 0 ) 

r' 


[ 1 ] 

-(T) 


+ O(^), 


(87) 


where Di is given by (^ and must satisfy the relation 


V. COMPARISON WITH FOUR-DIMENSIONAL EINSTEIN GRAVITY 

In this section we compare gravity in the brane world scenario with 4-dimensional Einstein gravity. We show that 
4-dimensional Einstein gravity is restored at low energies on our brane, provided that matter fields on the hidden 
brane are not excited. This qualitative result does not depend on which brane we are living. However, quantitatively, 
the effective 4-dimensional gravitational constants on the two branes are different. 


A. Perturbations in four-dimensional Einstein gravity 

For the purpose of comparison, we review linearized Einstein gravity in 4-dimensions in terms of gauge-invariant 
variables. We consider general perturbations around the Minkowski spacetime. Namely, we consider the metric 

dsl = (g^°J + 5q^^)dy^dy^, ( 88 ) 

where 

9 ® = (89) 

H is a non-zero constant conformal factor , and 

ddfiv = d(T)T{T)^iv + d(LT)^(LT)/ii/ + d + d(Y')T(Y)^iv (90) 

Here, the coefficients cf{t)i ’^{lt)i '^(v) are constants. Following the previous section, we can construct 

gauge-invariant variables /(t) and /. 

/(T) =d(T), 

f = a(y) -f ■ (91) 

As for the stress energy tensor we consider it as a first order quantity and expand it as follows. 

Sij.1, = -I- f(LT)T{LT)nv + ^LL)T(LL)iiv + d(jf)Ti^Y)fj.iy, (92) 

where the coefficients T(t lt,ll,y) are constants. In the Minkowski background these coefficients are gauge-invariant 
by themselves. 

We can expand the four-dimensional Einstein equation G(^Y}tJ.v = up to first order in the perturbations, 

where G( 4 )^,y is the Einstein tensor constructed from the metric and express it in terms of the above gauge- 
invariant variables. The Einstein equation is 
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( 93 ) 


2 t ( f ) = ?>'n^''k^KT(LL), 
f = -167rGArf^^f(iL) 


for scalar perturbations, 


''■(LT) — 0 (94) 

for vector perturbations, and 

= 16ttGnT(t) (95) 

for tensor perturbations. 

Note that in such a gauge 0 that d^^(5qay — qj^Jq^°'>^'^6qpcr/2) = 0, the variable / is given by / = (2/3)tT(y). Hence, 
in this gauge the second equation of (p3|) becomes 


q(°'>t^’'kpki,a^Y) = - IQttGn^y) 


(96) 


It is easy to confirm by noting □ = that equations (^) and (^) are equivalent to (|^) with B = 0. Of 

course it is also easy to derive Newton’s law from equation (3) with AK = B = 0 


B. Perturbations in the brane world scenario and restoration of Einstein gravity 


Here we shall apply the results of Sec. IV to the case that our 4-dimensional universe is one of the two branes S± 
and that there is no matter excitation on the other brane. 

First, let us suppose that our 4-dimensional universe is S+ and that there is no matter excitation on the other 
brane S_. For scalar perturbations, by setting T_(y) = 'j^-(LL) = 0, we obtain 


9 —24 I — 

- 4e 

J+ |■W+ ® 




(97) 


for zero modes. This equation can be considered the lowest order equation in an expansion in terms of /r = 
where I is a characteristic length scale of the model Moreover, we have the conservation equation 

2f+(y) = 34''kpk^f+,^LL) (98) 


for all modes. 

For vector perturbations, by setting = 0 we obtain 

^+{LT) = 0 (99) 

for general k^. Here, we have used the relation (^) for zero modes, or modes with q^'^k^ki, = 0. For tensor 
perturbations, by setting ~ obtain 


( 0 ) 


q+^'^^'kpk^fYiT) = 


2 — 2 ^ 4 +—..., 

nie + + 


dw' 


t+{t) 


( 100 ) 


for small but non-zero |/r| = 4ri^''kpk^\. Here, we have used the relation (81) to show that = 0. Actually, 


eq. (1^) remains correct even for rj^^^kpk, 


0 since D 2 in ( [751) is arbitrary and the relation & implies 'r_|_(y) = 0. 
Next, let us suppose that our 4-dimensional universe is S- and that there is no matter excitation on the other brane 
S+. Repeating the above analysis, we can obtain the same result with ± replaced by =F. For scalar perturbations, we 
have 


(3) and (4) we adopted this gauge choice. 

®See eq. (^^ for tensor perturbations and the footnote before it. 
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( 101 ) 


2T_(y) = 2,r]f^''kf,k^T_(^LL), 


I- 


^2-'2A-- 


.( 0 ) 


-2A--, 


.,( 0 ) 


dw' 

Jw_ 


T-(LL) 


for zero modes. 

For vector perturbations, we have 


'’'-(LT) = 0 


( 102 ) 


for general k^. For tensor perturbations, we have 


q-^^''k^kyf_{T) = 


2„-2A_-. 


do) 


Kte 


J w_ 


(103) 


for \fj,\ zero or sufficiently small. 

Therefore, by setting ^ 4-dimensional Einstein gravity is restored at low energies, provided that 

there are no tachyonic modes, or positive p, modes, with T±(r,= 0- The condition T±(r,LL,F) = 0 for scalar 
and tensor perturbations corresponds to the boundary condition 




] 


= 0 , 


w—w± 

{e^^F,T)y\ =0, 

'' \ W — WA- 


(104) 


where B± is given by Hence, we need to show the absence of tachyonic modes given this boundary condition. 

For this purpose we need to specify model potentials and background solutions since it seems difficult to obtain 
analytic solutions to the bulk equations (|^ and (^ for a general positive general potentials, and a general 

background. On the other hand, we have already shown that there are no tachyonic modes for vector perturbations 
(see eq. (^^). The effective 4-dimensional Newton’s constants Gn± on are given by 


167rGiv± 


Kte 


-2A±-af 


JW- 


(105) 


Since W- < w+, G±n is positive on both branes^ It is worth while mentioning that the following naive expectation 
holds. 


Gn+ 

Gn- 


„-2A+-, 


4°) 


„-2A_-a 


( 0 ) 


-(0) 
-( 0 ) • 


(106) 


where q‘£l^ 


are the unperturbed physical metrics on 




respectively. 


VI. SUMMARY 

We have investigated gravity in brane worlds. In particular, we considered a general model on a Z 2 -orbifold with 
a bulk scalar field and investigated gauge-invariant perturbations around a background with two parallel Minkowski 
branes. By using the doubly covariant formalism of perturbation of Israel’s junction conditions [Q developed in 
ref. p^ , we have obtained equations governing the low energy dynamics of doubly gauge invariant perturbations 
on the branes. These doubly gauge invariant perturbations include scalar, vector, and tensor perturbations. After 
that, we have shown that 4-dimensional Einstein gravity is restored on both branes at low energies, provided that 
combinations of the bulk and the brane potentials B± given by (|^) are not zero at the branes and that there are 


®If we had choose the convention W- > w+ then the sign in (|^) would be reversed and G±jv would still be positive, as should 
be. 
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no tachyonic modes for scalar perturbations with the boundary condition (104). It is easy to see that the condition 
B± ^ 0 is equivalent to the absence of a scalar zero mode with vanishing matter on the branes. Actually, the would-be 
zero mode with vanishing matter is prohibited by the boundary condition ( |60|) with t±(ll) = 0 if and only if B± = 0. 

The effec tive 4-dimensional Newton’s constants Gn± on S±, respectively, are given by (105) and the naive expec¬ 
tation (1^) holds. Thus we have confirmed and extended the result of that pure Einstein gravity is restored for 
the braneworld scenario with a single bulk scalar at low energies. 

The condition B± 7 ^ 0 is equivalent to the condition that the inter-brane distance (radion) is stabilized (to be 
precise, fixed) and that the background bulk scalar field has non-zero derivatives near the branes. To see this it is 
enough to differentiate the background scalar-field matching condition (the second of (^)) with respect to w± [p^ . 
In fact, the radion stabilization can be stated as 


0 7 ^ e 


^±. 


d 


dw± 


i\i>(o) 


1 


w—w± 








( 0 ) 


^,( 0 ) 




since the following equation is always satisfied because of the background equation. 


d 


dw± 


e^A 


1 


w—w± D 


= 0 . 


(107) 


(108) 


Essentially, the condition (IW) says that the position of the brane cannot be changed without changing the bulk 
solution. Hence, actually the condition B± 7 ^ 0 is equivalent to the condition that 7 ^ 0 and that the inter-brane 

distance is fixed. Since is related to by (p7| ) on branes, finally the condition B± 7 ^ 0 is equivalent to the 
condition that the background bulk scalar field is changing near the branes and that the radion is stabilized. 

One might expect that it is sufficient but not necessary to fix positions of both branes in order to fix the distance 
between two branes. In fact, it is necessary to fix positions of both branes as in (107). To see this, let us assume that 


the position of one brane can be changed without changing the bulk solution at all. (The position of the other brane 
can be either fixed or unfixed.) In this case it is impossible for the position of the unfixed brane to affect the position 
of the other brane since the only possible medium of the communication between the two branes would be the bulk 
but the bulk solution is independent of the position of the unfixed brane by assumption. Thus, the distance between 
the branes is not fixed in this case. In other words, the matching condition for the background with the 4-dimensional 
Poincare symmetry is local. 

Although to derive Newton’s gravitational couplings, we have used only matter source at the brane S-|_ (visible 
brane), it is instructive to write down the most generic weak gravity equations taking into account the contribution 
from the hidde n brane, which enters the final results in a symmetric manner, see the last equations in subsections 
Suppose energy-momentum tensor at the visible brane is and that on the hidden brane E_ is 


IV A and IV C 


Let us again, as in Sec. II, use linearized metric perturbations at the visible brane. Similar to (jg), in most 
general case we have 


ah = 167 rG +(1 - 2B+)T+ + 16nGjj{l - 2B-)T 


(109) 


Therefore, in principle, we need four gravitational couplings, two for tensor modes and two for scalar modes, to 
describe the linearized gravity at the visible brane. We also learned that each of them are constructed from the 
fundamental 5d gravitational constant, K 5 , and combinations of the warp factor and “effective” extra dimensional 
distance I ~ \ It turns out that the BD parameter, which characterizes the gravitational coupling for 

the scalar mode, depends on the bulk scalar sector. In terms of coefficient B of Sec. II, B^ is none-zero for the case 
without bulk scalar, B^=Q for the case of a single scalar field. We conjecture, based on our derivation with a single 
scalar, that in general case of more than one bulk scalar, B^ again will be non-zero. It will be interesting to calculate 
it explicitly for n ext simplest case of two bulk scalars. 

Equation (10£) may have interesting insights from point of view of superstring (Horava-Witten) phenomenology. 
Indeed, it says that there are, in general, four different Planckian masses! It also teaches us that derivation of 4d 
effective theory by means of f dw integration will wash out the subtleties of the actual gravitational dynamics. Finally, 
from point of view of gravity and its application to astronomy, it gives us rather unusual theory with four gravitational 
couplings. 
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APPENDIX A: HARMONICS IN MINKOWSKI SPACETIME 

In this appendix we give definitions of scalar, vector and tensor harmonics in an n-dimensional Minkowski spacetime. 
Throughout this appendix, n-dimensional coordinates are (/r = 0,1, • • •, n — 1), is the Minkowski metric, and 
all indices are raised and lowered by the Minkowski metric and its inverse . 


1. Scalar harmonics 


The scalar harmonics are given by 

Y = exp{—ikpx'’), (Al) 

by which any function / can be expanded as 

f = J dk cY, (A2) 

where c is a constant depending on k. Hereafter, k and dk are abbreviations of {k^} (/i = 0,1,- - - ,71 — 1) and 
respectively. We omit k in most cases. 


2. Vector harmonics 


In general, any vector field Vp can be decomposed as 

Vf, = V(T)p, + dpf, (A3) 

where / is a function and V(^t}ii is a transverse vector field: 

= 0. (A4) 

Thus, the vector held can be expanded by using the scalar harmonics Y and transverse vector harmonics V(t)/x 
as 


Vfj. = j dk [c(T) V(t)m + C(l)5^E] ■ 

Here, c^x) and C(i) are constants depending on k, and the transverse vector harmonics V(t)/j are given by 

ViT)!! = Ufj,exp{-ikpxP), 

where the constant vector Up, satisfies the following condition. 

k^Up = 0 


for k^kp ^ 0, and 


k^Up = 0 , 

T^Up = 0 


(AS) 

(A6) 

(A7) 


(AS) 


for non-vanishing kp satisfying k^kp = 0, where is an arbitrary constant timelike vector. For kp = 0, the constant 
vector does not need to satisfy any of the above conditions. For the special case k^kp = 0, the second condition in 
(AS) can be imposed by redefinition of C(i). Actually this condition is necessary to eliminate redundancy. Note that 
the number of independent vectors satisfying the above condition is n — 1 for k^kp ^ Q and n — 2 for k^kp = Q and 
that these numbers are equal to the numbers of physical degrees of freedom for massive and massless spin-1 fields in 
n-dimensions, respectively. 

Because of the expansion ([A5|), it is convenient to define longitudinal vector harmonics V(l)^ by 






= dpY = -ikpY. 


(A9) 
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3. Tensor harmonics 


In general, a symmetric second-rank tensor field can be decomposed as 

“b “b “b (AlO) 

where / is a function, is a vector field and is a transverse traceless symmetric tensor field: 


— 0 - (All) 

Thus, the tensor field can be expanded by using the scalar harmonics F, the vector harmonics V(^x) and 
and transverse traceless tensor harmonics T(^t) as 


= J dk -I- C(^LT)id^V(^T)u + duV(^T)fj.) 

+C(LL)(9^V(i)i, + -I- C(y)F77^;,] . 


(A12) 


Here, C( 7 ’), C(ll), and C(y) are constants depending on fc, and the transverse traceless tensor harmonics 

are given by 


T(^T)p.iy = Sfj.„exp{-ikpxf), (A13) 

where the constant symmetric second-rank tensor s^i, satisfies the following condition. 

k^Spi, = 0 , 

< = 0 (A14) 

for k^kp ^ 0, and 

k^Sfj,„ = 0 , 

50=0, 

= 0 (A15) 


for non-vanishing k^ satisfying k^k^ = 0, where is an arbitrary constant timelike vector. For k^ = 0, the constant 
ten sor S uu does not need to satisfy any of the above conditions. For the special case k^^kp = 0, the last condition 
in ( A15 ) can be imposed by redefinition of ^{ll) and C(y). Actually this condition is necessary to eliminate 

redundancy. Note that the number of independent symmetric second-rank tensors satisfying the above conditions is 
(n -I- l)(n — 2)/2 for k^kp ^ 0 and n{n — 3)/2 for k^k^ = 0 and that these numbers are equal to numbers of physical 
degrees of freedom for massi ve an d massless spin-2 fields in n-dimensions, respectively. 

Because of the expansion ( [A12| ), it is convenient to define tensor harmonics T'(lt)i and T(y) by 


T(LT)iiv = dpV(^T)i' + 

— %{TLpkjj Ui/kp^Y^ 

2 

T{ll)iju = dp,V(^L)u + di,V(^L)ti — 

= (^-2kpk„ + ^kPkpijpJj Y, 

T(Y)pv = ■npuY. (A16) 
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